Temperature dependence of the Kondo resonance and its satellites in CeCu 2 Si2 
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We present high-resolution photoemission spectroscopy studies on the Kondo resonance of the 
strongly-correlated Ce system CeCu2Si2. Exploiting the thermal broadening of the Fermi edge 
we analyze position, spectral weight, and temperature dependence of the low-energy 4f spectral 
features, whose major weight lies above the Fermi level Ef- We also present theoretical predictions 
based on the single-impurity Anderson model using an extended non-crossing approximation (NCA) , 
including all spin-orbit and crystal field splittings of the 4f states. The excellent agreement between 
theory and experiment provides strong evidence that the spectral properties of CeCu2Si2 can be 
described by single-impurity Kondo physics down to T ~ 5 K. 

PACS numbers: 71.27.+a 71.28.+d 79.60.-i 71.10.-w 



Many of the salient properties of intermetallic rare 
earth compounds originate from strong electronic correla- 
tions in the rare earth 4f orbitals. The hybridization with 
the conduction electron continuum leads to complex low- 
temperature behavior, summarized by the term Kondo 
scenario Q|. It includes quenching of the 4f magnetic 
moments below the Kondo temperature Tr, a strong en- 
hancement of thermodynamic quantities like the mag- 
netic susceptibility and the specific heat. The latter is 
ascribed to an enhanced electron density of states (DOS) 
at the Fermi energy Ef, the Kondo resonance (KR). In 
lattice systems, below a so-called lattice coherence tem- 
perature T* , the transport properties Q can often be de- 
scribed by the assumption of propagating quasiparticles 
with an effective mass enhancement of m*/m e = 10 2 to 
10 3 . The coherence temperature T* can be in the same 
range as Tk ||. Several of these heavy Fermion sys- 
tems Q undergo a magnetic ordering transition at low 
temperature T, become paramagnetic insulators (Kondo 
insulators) or even superconducting, like CeCu2Si2 [^J, 
where T K » 4.5 K to 10 K §, § and T* « 10 K §. In 
order to understand the interplay between local moment 
quenching and lattice coherence and ordering effects it is 
essential to know to what extent the system is described 
by local Kondo physics. However, it has been notoriously 
difficult to directly investigate the KR experimentally. In 
fact, the Kondo scenario as the origin of the 4f DOS en- 
hancement in Ce-based compounds has repeatedly been 
disputed in the literature (see || and references therein), 
mainly because the enhancement can still be seen at tem- 
peratures almost two orders of magnitude above Tr. 

There exist three established, spectroscopic meth- 
ods to investigate the spectral function close to the 
Fermi energy, photoemission (PES), inverse photoemis- 



sion (IPES) and scanning tunneling spectroscopy (STS). 
In the case of Ce compounds the 4f occupation is close 
to unity (4f : ), the Kondo temperature lies in the range 
Tr ps 1, . . . , 1000 K j|, and the KR has its maximum 
above Ep. The difficulty in observing the KR lies in the 
limited IPES energy resolution of O(100 meV), whereas 
PES has, in general, merely access to the tail of the KR 
below Ep and STS || probes the local 4f spectrum only 
indirectly by tunneling into the conduction states fief . 

In this Letter we present a careful analysis of high- 
resolution PES data close to the Fermi level, which 
allows to acquire information of the spectral function 
above the Fermi level. We restrict ourselves here to 
the heavy fermion compound CeCu2Si2, although exper- 
imental method and data analysis [ p~T| have also been 
successfully applied to other Ce intermetallics |Q. A 
comparison with calculations of the 4f spectral function, 
based on the single-impurity Anderson model (SIAM) 
including spin-orbit (SO) as well as crystal held (CF) 
splitting of the local orbitals, shows striking agreement 
between experiment and theory. In particular, we explain 
the temperature dependence of the spectral features up 
to T = 200 K. 

The PES experiments have been performed using a 
SCIENTA SES 200 analyzer and a monochromatized 
GAMMADATA VUV-lamp at photon energies of hv = 
21.2 eV (He I) and hi> = 40.8 eV (He II). A description 
of the experimental setup and the calibration of the ex- 
perimental parameters for the high-resolution measure- 
ments, like Ep, sample temperature T, and energy reso- 
lution AE, can be found in more detail in Ref. The 
poly-crystalline samples were cleaved in situ at low T to 
prepare clean surfaces. Because of the high surface re- 
activity of rare-earth compounds Jl4| the duration of a 
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FIG. 1: Theoretical 4f spectral function from calculations 
based on the SI AM using the NCA for T = 11 K and model 
parameters as in Fig. [| The hatched region is the occupied 
part of the spectrum. Energies are given in units of the half 
band width D — Ef- The inset shows the near-_EF region. 
The spectral features A - F are explained in the text. The two- 
electron state f 2 lies far outside the displayed energy range 
at « U + e/i. 
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measurement was kept less than about 12 h. 

In CeCu2Si2 the 7 spin-degenerate 4f levels are split 
by SO coupling into a total angular momentum J = 5/2 
sextet and an excited J = 7/2 octet. These are, in turn, 
CF split into 3 and 4 Kramers degenerate doublets, re- 
spectively. Thus, the single-impurity system is described 
by the multi-orbital Anderson Hamiltonian 



FIG. 2: Photoemission spectrum of CeCu2Si2 (hu — 40.8 eV, 
He II). Open circles represent the raw data close to the Fermi 
level, including the tail of the KR and the low-energy CF exci- 
tations at approximately 21 meV (AE = 6 meV, T = 11 K). 
Filled circles represent the raw data divided by the experi- 
mentally broadened FDD (see text). The shaded area marks 
the unreliable spectral range above 5ksT. The inset shows 
an extended energy range (measured with AE ~ 13 meV, 
T = 17 K) including the spin-orbit (J=7/2) excitation at 
« 270 meV below Ef and the non-4f background as deter- 
mined using He I (hu = 21.2 eV) radiation. 
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where Hn = Y% erd c- describes the conduction band 



with dispersion and creation operators ct for elec- 
trons with spin a. Sf m < Ep, m = 1,...,7, are the 
SO and CF split 4f single-particle levels with the corre- 
sponding creation operators f^ na . The hybridization ma- 
trix elements lead to an effective coupling matrix be- 
tween 4f states, T mn = tt V m% A k^ V kn^ where A k^°) 
is the conduction electron spectral function. The local 
Coulomb repulsion, known from IPES 15l 16) , is substan- 
tially larger than and, hence, effectively suppresses 
any double occupancy of the Ce 4f levels (U — * oo). 
In the multi-orbital system the Kondo temperature is 
given by T K « V2JE7 exp[-l/(2iV(0) J)], where N(0) 
is the conduction electron DOS at the Fermi level and 
the Kondo coupling J is obtained from a Schrieffer- Wolff 



transformation including the CF and SO states as 
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It is seen that J, and hence Tk, is enhanced due to cou- 
pling to the CF and SO split 4f states. 

To achieve a detailed comparison with experimental 
results, we have calculated the 4f spectral function of 
the model Eq. (|l|) using the non-crossing approximation 
(NCA) 111 including all CF and SO excitations. The 



NCA is known to reliably describe the low-temperature 
scale Tk and the position, the spectral weight and the 
life-time broadening of the peaks in the local spectral 
density down to T w 0.1 Tk 0- The spectrum has 
generically six distinct features as shown in Fig. [l| (A-F). 
They can be understood as follows: At low T the occupa- 
tion of the lowest 4f level is close to unity (n/i < 1), while 
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all other single-particle 4f states are essentially empty 
(n/ m w 0, m = 2, . . . , 7) because of the strong Coulomb 
repulsion £/. Hence, the broad 4f 1 — * 4f° ionization 
peak (A) with a full width at half maximum (FWHM) 
of r w Ti m corresponds to the lowest single-particle 
level £fi. Resonant spin flip scattering of electrons at the 
Fermi energy induces the narrow Kondo resonance (D) of 
width ~ ksTK, somewhat shifted above Ef due to level 
repulsion from the single-particle levels e/ m < Ef- The 
SO and the CF satellite peaks appear in pairs (B, F) 
and (C, E), respectively. They arise from virtual tran- 
sitions from the ground state into the (empty) excited 
SO (F) and CF (E) states and vice versa (B and C). The 
positions of the satellite peak pairs are, therefore, ap- 
proximately symmetrical about Ep. However, while the 
features above Ep have significant spectral weight, those 
below Ep appear merely as weak shoulders. This is be- 
cause the transition probabilities carry a detailed balance 
factor w — n l (l — n?), where ri 1 (n/) is the occupation 
number of the 4f orbital in the initial (final) state, i.e. w 
is large for the excitations E, F, but small for the transi- 
tions B, C. 

As n l , w are controlled both by T and U, the CF and 
SO satellites are signatures of strong correlations and are 
T-dependent. All these features are well described by the 
NCA (Fig. |l|). We note in passing that similar spectra 
should be observed in multilevel quantum dots in the 
Kondo regime. 

In order to identify the low-energy features in ex- 
perimental spectra as due to the physics of the single- 
impurity Anderson model (Eq. [l]), it is highly desirable 
to gain access to the sharp resonances D, E above Ep. A 
typical high-resolution spectrum on CeCu 2 Si2 is shown 
in Fig. ||, measured at T — 11 K with an energy resolu- 
tion of AE = 6 meV. The ionization peak A (not shown) 
appears at a binding energy of « 2 eV below Ep, the 
SO satellite B at « 270 meV. Fig. | shows the near- Ep 
region with the tail of the KR and a distinct spectral 
feature at about 21 meV, which can be assigned to the 
CF splitting of the J = 5/2 4f levels @. In this plot 
the spectral intensity at energies of a few meV above 
Ep and higher is completely suppressed by the Fermi- 
Dirac distribution (FDD). IPES experiments performed 
on similar systems do not resolve the KR and its satellite 
excitations expected above Ep, but show only one broad 
structure near Ep and a second one at much higher en- 
ergies, corresponding to the doubly occupied IPES final 
state 4f* -> if [^M- 

Greber et al. p0[ have demonstrated that a careful 
analysis of PES data based on a division by the respec- 
tive FDD allows to investigate the spectral function up 
to energies of approximately 5fcsT above Ep, provided 
that the noise of the data and the experimental broaden- 
ing are small enough. We applied this procedure to the 
spectrum in Fig. ^| where we divided the raw spectrum 
by the FDD at T = 11 K convoluted with the spectrome- 
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FIG. 3: Upper panel: Theoretical T dependence of the 4f 
spectral function of CeCu2Si2. The inset shows the calcu- 
lated spectrum below Tk (see text) . Used model parameters: 
e/i = -2.4 eV, E F = 4 eV, CF splittings of the J = 5/2 sextet 
Acf = 30 meV and 36 meV, SO splitting A S o = 270 meV, 
hybridization V = 200 meV. Lower panel: Photoemission 
spectra for T = 15 K, 47 K, 105 K, and 200 K. The experi- 
mental spectra, divided by the FDD, are clipped on the high 
energy side at m 5/dbT. The inset shows the data on the same 
energy scale prior to the division. All spectra are normalized 
to the same intensity at w 100 meV and offset for clarity. 



ter function, described by a Gaussian of full width at half 
maximum (FWHM) AE — 6 meV. Obviously there ap- 
pears a narrow peak with a FWHM of 6 meV and a maxi- 
mum at about 3 meV above Ep- This is sufficiently below 
5ksT w 7 meV at which energy the scatter of the data 
gets large due to the exponential increase of the inverse 
FDD. The spectrum of the CF satellite at w 21 meV be- 
low Ef remains unchanged by this procedure. It should 
be mentioned that the described procedure is not exact 
and therefore can lead to artifacts (of the order of AE) if 
the temperature broadening is comparable with the en- 
ergy resolution. However, in the present case this is only 
relevant for the 11 K-spectrum which appears slightly 
shifted away from Ep. 

In the following we analyze the T dependence of the 
near-_Ep spectral features. Fig. pi lower panel, shows 
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the experimental data at several temperatures, divided 
by the broadened FDD as described above. The up- 
per panel of Fig. |^ displays the NCA spectral func- 
tions for the same temperatures. The model parameters, 
i.e. the single-particle energies e/ m and the hybridiza- 
tion strengths T mm i, have been taken from independent 
experimental measurements where possible [ fl9| , |22] ] and 
have only slightly been adjusted (see caption of Fig. ||) to 
get optimal quantitative agreement with our experimen- 
tal results. The experimentally observed temperature de- 
pendence of the spectra is accurately described by the 
theoretical simulations, and a comparison allows to iden- 
tify the origin of the various spectral features: As already 
seen in Fig. ||at low T (T = 15 K) the KR at E F appears 
as a narrow line with a FWHM of about 6 meV. Due to 
the given energy resolution of AE = 6 meV the intrinsic 
linewidth of the Kondo resonance must be much smaller 
than AE at low temperatures. In order to estimate Tk 
from the data we have determined the model parameters 



physics, although the Kondo local moment quenching 
and lattice coherence effects occur on similar energy 
scales (T* Tk)- In particular, the Kondo resonance 
could be unambiguously identified, with a Kondo tem- 
perature in reasonable agreement with estimates from 
thermodynamic measurements. The persistence of an 
enhanced density of states even at room temperature is 
explained by the presence of crystal field excited states 
close to the Fermi energy. 
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to fi t the experimental results at their respective T and 
the n calculated the spectrum at T — 0.1 Tk ~ 1 K (inset 
of Fig. ^, upper panel), i.e. close to the unitarity limit 
which is not directly reachable in the experiment. Tk was 
then determined from the peak width of this theoretical 
low-T spectrum as Tk ~ 6 K. This result is in remark- 
able agreement with thermodynamic bulk measurements 
H 0, M , considering especially that the spectral function 
— and therewith Tk — as measured by surface sensitive 
PES might be modified by surface effects, as pointed out 
by several authors [|l4|, ^l], ^| . Towards higher T, the line 
width of the KR increases, while the maximum intensity 
becomes smaller and the distinction of KR and CF satel- 
lite is successively smeared out. At T = 200 K and above 
the thermal broadening of the FDD is large enough for 
PES to have access to the CF excitations at ~ 40 meV. 
At these elevated T the KR has disappeared as a sepa- 
rate peak, and the enhanced 4f DOS near Ep is rather 
due to the CF excitations which are broadened and posi- 
tioned somewhat above Ep, as seen in Fig. |. Thus, the 
persistance of an enhanced DOS at the Fermi level even 
at room temperature, despite a substantially lower Tk, 
is naturally explained within the SIAM in combination 
with CF excited states. 

In conclusion we have demonstrated that high resolu- 
tion photoemission spectroscopy gives detailed access to 
the spectral function of heavy Fermion Ce compounds 
and allows to investigate the temperature dependence 
of the Kondo resonance. Within an energy range of 
5ksT above the Fermi level several structures can be 
resolved, whose origin and temperature dependence are 
consistently described within the single-impurity Ander- 
son model. This excellent agreement provides strong ev- 
idence that the 4f photoemission spectra of CeCu2Si2 
are an immediate consequence of single-impurity Kondo 
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